Phase-coherent effects in multiterminal superconductor/normal metal mesoscopic 

structures 



O 

o 
o 



(N 



S 

T3 



- 1—^ 

X 



A.F. Volkov* 1 * and R. Seviour*, 
* School of Physics and Chemistry, 
Lancaster University, Lancaster LAI J^YB, U.K. 
^Institute of Radioengineering and Electronics of the Russian 



■ Academy of Sciencies, Mokhovaya str.ll, Moscow 103907, Russia. 

(February 1, 2008) 

In this report we analyse three effects which may arise in a mesoscopic multiterminal S/N 
structure (two normal and two superconducting reservoirs). We show that the Josephson critical 
current is a non-monotonic function of a voltage between the normal reservoirs. The influence of the 
spin-polarized electrons on the Josephson critical current was also studied. We show that if there 
is a temperature difference between the normal reservoirs, a voltage between the superconductors 
and normal conductors arises which oscillates as the phase difference varies. Its magnitude is much 

■ greater than the thermoemf in the case of the ordinary thermoelectric effect. 

o 

. 1. Introduction 

In this report we present results of the theoretical study of transport in diffusive 4-terminal superconductor-normal 
metal (S/N) mesoscopic structures. A great deal of attention was paid recently to the problem of the influence a 
current between the N reservoirs has on the Josephson effect in the superconducting curcuit. This problem was 



analysed first in Ref. JjJ for a ballistic S/N system and in Ref. ]2| for a diffusive system. It was shown that the 
Josephson critical current I c changes sign at a certain voltage V between the N reservoirs, and the Josephson junction 
switches into a 7r-state (7r-junction). The reason for the sign reversal effect is that the distribution function in the N 
wire, which can have an equilibrium form, is shifted with respect to the distribution function in the superconductors 
by the value V. Later the dependence I C (V) was calculated in various models and approximations in Refs. 
(diffusive case) and P,|lOt (ballistic case). The sign reversal effect in a diffusive Au/Nb mesoscopic structure was 
observed by Baselmans et al. |TJJj . 

In multi-terminal S/N structures one can observe not only the sign reversal effect, but also a number of other 
interesting phenomena. It was shown in Refs. || that the measured critical current I m in a structure similar to 
that shown in Fig. la may not only decrease, but also may increase with increasing voltage V . In particular one can 
observe the Josephson-like effects (plateau on the Iz{Vs) curve, oscillations of the measured critical current I m in a 
' magnetic field etc) even if the Josephson coupling between superconductors under equilibrium conditions is vanishing. 
£C) . The stimulation of the Josephson effect in this structure is related in particular to a phase dependence of the S/N 
interface or N wire conductance: 5G ~ cos tp. It turns out that apart from the Josephson term I c sin tp, an additional 
term I sg cos tp appears in a dynamic equation for tp. This term leads to a change of the critical current from the value 

I c to a new (measured) value I m = J 1% + I% g , where I sg is the amplitude of the subgap current. The stimulation of 

the Josephson effect by a current between the N reservoirs was analysed in Refs. || on a simplified model of gapless 
{""] superconductors. Here we present results of the analysis for the real case where the superconductors have a gap and 
the density-of-states (DOS) has a singularity Qj. We show that the measured critical current I m reaches a maximum 
at eV = A and depends weakly on the temperature T. These results agree with recent observations of the stimulation 
Q | of the Josephson critical current E^] . 

O ■ We also study the effect of the spin injection into the N wire from the ferromagnetic (F) reservoirs on the critical 
\ current I c for the structure shown in Fig. lb. If the magnetizations in the F reservoirs have the same direction, then 
there is no spin injection. In the case of antiparallel orientation, the total magnetization of the injected electrons in 
the N wire is finite. We show that in the case of a N wire short length 2L (L < {L e , L sp }; here L e and L sp are the 
■ energy and spin relaxation length) the spin injection does not affect the critical current and the dependence I C (Y) 
is the same as in the case of normal (nonmagnetic) reservoirs. In a semi-mesoscopic structure (L e < L < L sp ) the 
7r-junction is realised only in the case of the antiparallel orientation of the magnetizations in the F reservoirs when 
spins are accumulated in the N wire. 

Finally we analyse the thermoelectric effect in the structure shown in Fig. la. We assume that there is no current 
between the N reservoirs, but the reservoirs are maintained at different temperatures T = T Q ± 5T. It will be shown 
that the temperature gradient leads to a voltage difference Vr between the normal and superconducting circuits. This 
voltage does not contain the small parameter Tjtp (as is the case for the ordinary thermoeffect in normal metals) 
and reaches a value of order 5T/e. In addition, the voltage Vt depends on the phase difference tp; it is zero at tp — 
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and oscillates with increasing ip. 
2. Basic equations 

The Green's function technique is perhaps the most convenient and powerful method for studing transport and 
nonequlibrium phenomena in S/N mesoscopic structures. If we are not interested in purely quantum interference 
effects arising as the result of multiple (non-Andreev) reflections, then we can use a simpler, well developed method of 
quasiclassical Green's functions (see, for example, ]i"3| ). In this method the matrix Keldysh function G is introduced 
along with the retarded (advanced) matrix Green's functions G R ( A > = G r ^t z + F r ^ a \ here F R ^ is the matrix 
condensate Green's function. The functions G^" 4 ' determine the excitation spectrum and the DOS. The function G 
describes the nonequilibrium properties and is related to the matrix distribution functions / : G a — G R f a — f a G A , 
where f a = / Q +l + fa-%- The components of / a can be expressed through the distribution functions of electrons 
n a and holes p a (a is the spin index). The function f a+ describes an electron and hole distribution symmetrical in 
branch populations: f a+ = 1 — (n a +p^), here a =J, if a =|. The function / a _ = — (n a — p&) describes the branch 
imbalance and determines the electric potential. 

The equations for G a and Gf' A ' can be solved analitically in the case of a short structure (cth = D/L 2 >> e, 
where e is a characteristic energy; e = min{T, A}) p4|-p8[ . In the general case one needs to solve equations for the 
distribution functions f a ± and the Usadel equation |17| , |l9| -|2l[|. In the present paper we consider the case when the 
temperature may be both greater or less than the Thouless energy exh and the approximation of a short length is not 
valid. The distrubution functions in the structure shown in Fig. la obey the kinetic equation (we consider only the 
dirty limit) 

Ld x [M±d x f±(x)+JsMx) ± Jand x f T (x)}= r[ A± 5{x - L x )+ A± S(x + LJ]. (1) 
where M± = (1 - G R G A T (F it F A ) 1 )/2; J an = {F R F A ) z /2, J s = {l/2){F R d x F R - F A d x F A ) Zl A± = (vv s + 

sit)(/± -fs±) - (g z ±fs T +g ZT f T );gi± = (i/*)[(F R ±F A )(F R ±F A )] 1 ;g z ± = (i/a)[(f r t f a )(f r ±f a )] z ; 

The parameter r = R/Rb is the ratio of the resistance of the N wire R and S/N interface resistance Rb', the 
functions A- and A+ coincide with A-,A+ if we make a substitution tp — > —ip. We introduced above the following 
notations (F R F A )! = Tr(F R F A ), (F R F A ) Z = Tr(? z F R F A ) etc.; u, v s are the DOS in the N film at x = Li and in 
the superconductors. The functions fs± are the distribution functions in the superconductors which are assumed to 
have equilibrium forms. This means that fs+ = f e q = tanh(e/3) (here (3 = 1/(2T)) and fs- — 0, because we set the 
potential of the superconductors equal to zero (no branch imbalance in the superconductors). We assumed that the 
width of the S/N interface w is much less than and introduced the 5— functions. Note that in some papers the 
last term in the left-hand side of Eq.(l) is missing. One can solve Eq.(l) and express the distribution functions f a ± 
in terms of the retarded (advanced) Green's functions which should be found from the Usadel equation. Integrating 
Eq.(l) once, we obtain 

M±d x f±(x)+Jsh(x) ± J an d x f T (x) = J 1)2± (2) 

where the indeces 1,2 relates to the intervals (0, L\),{L-y, L) and the constants Ji^± are partial flows. For example, 
the partial flow Ji,2- determines the electrical current in the intervals (0, L±),(L±,L) 

/■oo 

h, 2 = o I deJi >2 - (3) 
Jo 

Here a is the conductivity of the N wire. As it follows from Eq.(l), the flows Ji,2± are connected with each other 

J2± - Ji± = rA± (4) 

The functions rA± are the flows through the S/N interface. In particular, rA_ is the partial electrical current 
through the S/N interface. The first term (vvs + <7i+)/- m ^ ne expression for A- is the quasiparticlc current (the 
term g\+f- contributes to the subgap current), and the second term (g z +fs+ + 9z—f+) is related to the supercurrent. 
One can show that the partial supercurrent J$ is a constant in the interval (0, L\) and equals zero outside this interval. 

We can solve Eqs.(2) and express the distribution functions f± through the constants J\^±- At x — ±L the 
functions f± are equal to equilibrium distribution functions in the N reservoirs. Therefore we can find the functions 
f±(x) and substituting them into Eq.(4) determine f±(Li) and J\^±- Thus the problem is reduced to solving the 
Usadel equation. The solution of the Usadel equation can be found either numerically or analitically in limiting cases. 
In the next sections we will employ this method to study the dependence of the critical current on the voltage between 
the N reservoirs and the thermoeffect. 
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3. Suppression and enhancement of the critical current 

In this section we consider a symmetrical structure shown in Fig. la. We assume that the electric potential at 
the N reservoirs is V and at the superconductors is zero (the quasiparticle currents flow from the N reservoirs to 
the superconductors and the supercurrent flows between the superconductors). The distribution functions in the 
reservoirs have the equilibrium form: Fy± = [tanh((e + eV)(3) ± tanh((e — eV){3)]/2. 

First we consider the case of a small r, that is, the interface resistance is large and the condition 

(e^)- 1 << r « 1 (5) 

should be satisfied. The first condition ensures that inelastic collision term can be neglected (here r c is the energy 
relaxation time, e ~ min{T, tTh})- Then the distribution function / + is constant and in the main approximation is 
/+ = (F\/ + + f eq (r2W s )) / + i r 2vv s ) 1 here = rL^jL. The function /_ is constant in the interval (0, L\) and in 
the main approximation equals /_ = Fy— /(l -\-Tiw^). Outside this interval the function /_ increases linearly to the 
value Fv— • Following the method presented above, we find the current through the S/N interface 

h {V) = h - h = Io(V) + I sg {V)cos<p-I c {V) sin ip (6) 

The first two terms are the quasiparticle current and the last term is the Josephson supercurrent. All the components 
of the current depend on V. This expression shows that at a given control voltage V and zero voltage difference between 
superconductors (ip is constant in time) the current ^3 may vary with changing (p in some limits: | I%{V) — I (V) |< 
I m (V). This means a plateau on the Vsih) characteristics (see ||); here Vs — (h/2e)d t <p is the voltage difference 
between superconductors. We can write the phase-dependent part of I3 in the form — I m sm(ip + a), where 

I m = \J Ic + 1% IS the measured critical current, cos a = I c /I m . In Fig. 2 we plot the dependence of the measured 
critical current I m on the control voltage V for different temperatures. It is seen that the temperature dependence 
of the maximum of I m which is achieved at eV ~ A is much weaker than the I C (T) dependence. These results 
qualitatively agree with the experimental data |12| . 

4. Spin injection and the critical Josephson current. 

In this section we present results of the study of spin injection on the critical Josephson current I c for the structure 
shown in Fig. lb p2| . We calculate the distribution function / + which determines I c for parallel (ft) and antiparallel 
(11) orientations of the magnetization in the ferromagnetic (F) reservoirs. It is assumed that the voltage 2V is applied 
between the F reservoirs and a current flows between the N reservoirs. We consider different limits of the length of 
the N wire. 

4-.1. Mesoscopic limit: L < {L e ,L sp }. 

In this limit the function /+ = (f-\+ + fi+)/2 does not depend on mutual orientation of ferromagnets and has the 
same form as in the nonmagnetic case: /+ = Fy+. The conductance G and magnetic moment of injected spins M (per 
unit volume of the N wire) depend on mutual orientation of the ferromagnetic domains. In the case of the parallel 
orientation (ft) we obtain G = (2R^ + Rl)^ 1 + {2Rb[ + Rl)^ 1 and M — 0. In the case of the antiparallel orientation 
(U) we have G = 2/{R b] + R H + R L ) and M = 4 M (i? bT - R H )VP/a 2 (2R M + R L )(2R H + R L ){R b] + R 2l + R L ). 
Here i?*T,J. and Rl are the F/N interface and the N wire resistances for electrons with up and down spins . Since 
the interface conductances R^ 1 ^ differ from each other (they are proportional to different density of states in the 
ferromagnets Vp^^), spin injection takes place in the antiparallel (fj) configuration. However the critical current I c is 
the same for both configurations. Its dependence on V is similar to that found for the structure under consideration 
with normal (nonmagnetic) reservoirs ||. 

4-1. Semi-mesoscopic limit: L e < L < L sp . 

In this limit there is no spin relaxation in the N wire and the magnetization M is determined again by formulae 
presented above. On the other hand the magnetization is related to a potential V sp which determines an imbalance 
between the spin subsystems in the N wire: M = 2fiveV sp . The function / + = (/f+ + f[+)/2 has an equilibrium form 
corresponding to the potential V sp : /+ = [tanh(e + eV sp )(3 + tanh(e — eV sp ){3}/2. Therefore in the case of parallel (ff ) 
orientation I c does not depend on V as V sp = (V is proportional to V sp ). In the case of antiparallel (f j) orientation 
I c does depend on V sp in the same way as in the section 4.1. A similar effect of spin injection into a supereconductor 
on the energy gap A was analysed in Ref . p3| . The critical Josephson current is more sensitive to the form of the 
distribution function than the energy gap A. Therefore the effect considered here allows one to realise the it— junction 
and to make certain conclusions about the energy and spin relaxation rate. 

5. Giant thermoelectric effect 

It is well known that if the terminals of a normal conductor are maintained at different temperatures ST = T a ± ST, 
then in the absence of the current an thermoemf V em f appears at the terminals. The magnitude of V em f is equal 
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to: eV em f — Ci(T/cf)5T, where c\ is a factor of the order 1 and is the Fermi energy. In this report we study 
the thermoelectric effect in the mesoscopic structure shown in Fig. la. It will be shown that a temperature difference 
between the normal (N) reservoirs leads to a voltage between normal and superconducting curcuits Vr- The magnitude 
of this voltage does not contain the small parameter (T/cf), and besides it oscillates with a variation of the phase 
difference ip between the superconductors. 

We assume that the superconductors are connected via a superconducting loop and the phase difference between 
them ip is controlled by an applied magncitic field. The N reservoirs are disconnected and maintained at different 
temperatures T(±L) = T Q ± ST. We will calculate the electric potential in the N film and in particular the potential 
Vr in the N reservoirs . Since we set the potential in the superconductors equal to zero, the potential Vr is the voltage 
difference between the N reservoirs and superconductors which arises in the presence of the temperature difference 
ST. In the limit of high interface resistance (see condition (5)) we obtain from Eq.(2) in the main approximation in 
r : f+(x) = J + x + /+(0), where /+(0) = f eq = tanh(/3 e), and J+ = 5/3e/L cosh 2 (/3 D e), 6/3 = -/3ST/T . We assumed 
the ratio ST/T D to be small and expant the distribution functions in the N reservoirs. 

In order to find the function /_, we consider Eq.(2) and take into account that J2- is zero (no current through the 
N reservoirs). In the absence of the temperature gradient we obtain from this equation 

/_ = 0, Jsfeq = Jleq = r(g z - + 9z+)feq (7) 

The temperature gradient leads to the non-zero function /_ (x) which determines the electrical potential Vr (x) . We 
find from Eq.(2) 



f-( X ) = /_(0) + SJ lX + f dXiiJandJf+iX!) - J,Sf + ( Xl )] 

Jo 



(8) 



where 5f+(x) = f+(x) — f eq . From Eq.(l) we obtain for 5J\ = J\ — J\ eq 

SJ 1 = ±r[ff, + J/+(±Li) - 0i+/_(±Li)] (9) 

Here the signs ± relate to the points ±L\. From Eqs.(8) and (9) it follows that SJ\ = 0, and in the main approximation 
in r the distribution function is /-(0) « 8f + (±L 1 )g z+ / g 1+1 where 8f + {±L\) = ±J + L\. The electric potential Vr is 
related to the function /_ and is an even function of x 



L18T a f°° ^ g z+ e 

9i+ cosh 2 (/3 e) 



(10) 



Here we substitute the function /_(0) which is the main contribution to the function f-(x). Therefore the potential 
Vr arising in the presence of the temperature gradient is almost constant along the N wire. This potential equals 
approximately the voltage difference Vr between the N reservoirs and superconducting loop. It is worth noting that 
Vr determined by Eq.(10) does not depend on the small parameter r (however one should have in mind that the 
condition (5) imposes limits on this parameter). We calculated the integrand using both analytical an numerical 
solutions of the Usadel equation. In Fig. 3 we show the temperature dependence of Vr for = ir/2. We see that the 
dependence Vr(T) is nonmonotonic (reentrant behaviour). One can easily estimate Vr on the order of magnitude. 
We obtain 

eV T = oT^/L) sm <f i ) . ' 2n x ,\ (11) 

Here Ci^{f) are periodic functions of the phase difference ip; they are of order 1 and are not zero at <p = 0. If the 
ratio T/erh is of order 1, then the thermocmf is of the order ST, that is, several orders greater than the thermoemf 
in the normal metals. 

The physical explanation of the effect is the following. The temperature gradient creates a deviation of the distrub- 
tion function Sf+ = —S(n + p) from equilibrium. On the other hand the superconductors do not affect this function 
because complete Andreev reflections conserve the total number of excess electrons and holes. The function 6f+ has 
different signs at ±Li and leads to an additional Joscphson current SJs = rg z+ Sf + (Li) of the same sign at ±Li(sign 
of the function g z+ is different at these points). In order to compensate these additional currents, the potential Vr 
arises in the N wire producing a subgap current SJ sg — rg\ + f-(±L\) which compensates the current 5Js- 

We have neglected the ordinary thermoelectric current It in the N wire because it contains a small parameter (T/er) 
and leads to a small contribution to Vr - Recently the influence of the proximity effect on the ordinary thermoelectric 
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effect in the mesoscopic S /N structure was studied both theoretically |^4j and experimentally (25) . The effect analysed 
here is much stronger. 

In simplified models we have analysed three possible effects in the 4-terminals S/N mesoscopic structures; one of 
them (see section 3) has been observed experimentally. Further theoretical studies of these and other effects are 
needed to account for processes which were neglected in our analysis. It is interesting, for example, to investigate how 
the obtained results are changed by the energy relaxation in the N wire which seems to be faster than it was expected 

H ' 
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FIG. |l| Schematic view of the 4-terminal S/N/S structure under consideration. 




FIG. g The measured critical current (I m ) vs V for different temperatures: (3 = erh/2T. The parameters are: A 
10e Th ,L 1 /L = 0.3, r = 0.3. 




FIG. g The dependence of the normalised voltage Vr on inverse temperature j3 at ip = it/2 (the parameters are A/exh 
W,Li/L = 0.5). 
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